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PREFACE 


This  Report  is  one  of  a  series  describing  the  results  of  numerical 
experiments  with  the  two-level  Mintz-Arakawa  atmospheric  model.  These 
experiments  are  being  carried  out  as  part  of  the  Rand/ARPA  Climate  Dy¬ 
namics  Program,  whose  primary  concern  is  the  analysis  and  validation  of 
the  circulation  changes  resulting  from  selected  alterations  of  the  model 
atmosphere's  initial  and/or  boundary  conditions.  It  is  also  necessary 
to  compare  the  climatic  simulations  with  observations,  and  the  present 
work  is  an  initial  contribution  in  this  direction. 

Meteorological  studies  suggest  that  technologically  feasible  opera¬ 
tions  might  trigger  substantial  changes  in  the  climate  over  broad  re¬ 
gions  of  the  globe.  Depending  on  their  character,  location,  and  scale, 
these  changes  might  be  both  deleterious  and  irreversible.  If  a  foreign 
power  were  to  bring  about  such  perturbations  either  overtly  or  covertly , 
either  maliciously  or  heedlessly,  the  results  might  be  seriously  detri¬ 
mental  to  the  security  and  welfare  of  this  country.  So  that  the  United 
States  nay  react  rationally  and  effectively  to  any  such  actions,  it  is 
essential  that  we  have  the  capability  to:  (1)  evaluate  all  consequences 
of  a  variety  of  possible  actions  that  might  modify  the  climate,  (2)  de¬ 
tect  trends  in  the  global  circulation  that  presage  changes  in  the  cli¬ 
mate,  either  natural  or  artificial,  and  (3)  determine,  if  possible, 
means  to  counter  potentially  deleterious  climatic  changes.  Our  pos¬ 
session  of  this  capability  would  make  incautious  experimentation  un¬ 
necessary,  and  would  tend  to  deter  malicious  manipulation.  To  this  end, 
the  Advanced  Research  Projects  Agency  initiated  a  study  of  the  dynamics 
of  climate  to  evaluate  the  effect  on  climate  of  environmental  pertur¬ 
bations.  The  present  Report  is  a  technical  contribution  to  this  larger 
study. 

Other  Rand  publications  related  to  the  present  Report  include 
R-877-ARPA,  which  presents  the  detailed  documentation  of  the  atmo¬ 
spheric  model  employed. 
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SUMMARY 

\ 

The  global  distributions  of  the  mean  January  surface  wind  stress, 
the  net  di  tbatic  heating  rate,  and  the  net  rate  of  moisture  addition 
as  simulated  in  a  30-day  integration  with  the  two-level  Mintz-Arakawa 
atmospheric  general  circulation  model  are  presented.  The  latitudinal 
distributions  of  the  zonal  averages  of  these  forcing  fields  are  shown 
to  be  in  reasonable  agreement  with  the  available  observations.  The 
most  prominent  discrepancies  are  evidently  due  to  the  model's  simula¬ 
tion  of  excessive  convective  precipitation  (and  the  associated  con¬ 
vective  latent  heating)  in  the  tropics,  especially  in  the  northern 
(winter)  hemisphere.  The  zone  of  simulated  tropical  precipitation 
extends  some  15°  poleward  of  the  observed  position,  and  results  in  a 
corresponding  distortion  of  the  field  of  evaporation  minus  precipita¬ 
tion  (or  moisture-addition  rate). 

In  determining  the  monthly  mean  forcing  fields,  the  (convective) 
precipitation  rate  in  particular  should  be  accurately  accumulated  dur¬ 
ing  the  integration;  the  customary  use  of  6-hour ly  fields  results  in 
a  sampling  error  as  large  as  25  percent  for  even  the  zonally  averaged 
rainfall.  With  the  exception  of  a  small  sampling  error  in  the  mean 
rate  of  absorption  of  solar  radiation,  the  components  of  the  other 
forcing  fields  are  satisfactorily  determined  by  6-hourly  data. 


Preceding  page  blank 


-vii- 


ACKNQWLEDGMENTS 


Thanks  are  due  Professors  Yale  Mintz  and  Akio  Arakawa  of  the  De¬ 
partment  of  Meteorology,  University  of  California  at  Los  Angeles,  for 
permission  to  use  their  two-level  general  circulation  model  for  the 
present  simulation.  I  would  also  like  to  thank  Professor  Mintz  for 
his  constructive  review  of  the  present  paper.  Acknowledgment  is  due 
Robert  Mobley  and  Al  Nelson  of  Rand  for  their  supervision  of  the  com¬ 
putations  and  machine  graphical  output. 


Preceding  page  blank 


-ix- 


CONTENTS 

PREFACE . iii 

SUMMARY  .  v 

ACKNOWLEDGMENTS . . 

Section 

1.  INTRODUCTION  .  1 

2.  THE  FORCING  FIELDS  IN  THE  MODEL .  3 

3.  THE  GLOBAL  DISTRIBUTION  OF  THE  SIMULATED 

FORCING  FIELDS . 5 

A.  The  Net  Frictional  Force,  ? .  5 

B.  The  Net  Diabatic  Heating  Rate,  H . 10 

C.  The  Net  Rate  of  Moisture  Addition,  E  -  P . 13 

4.  THE  ZONAL-AVERAGE  SOURCE  TERMS  AND  THEIR 

COMPARISON  WITH  OBSERVATION  .  17 

A.  The  Surface  Frictional  Stress  Components  .  17 

B.  The  Components  of  the  Net  Heating  Rate . 19 

C.  The  Components  of  the  Moisture  Balance . 26 

5.  SUMMARY  AND  CONCLUSIONS . 32 

Appendix 

THE  EFFECTS  OF  SAMPLING  FREQUENCY  ON  THE  ESTIMATION 

OF  THE  MEAN  FORCING  FIELDS . 34 

REFERENCES . 43 


Preceding  page  blank 


-1- 


1.  INTRODUCTION 


The  two-level  Mintz-Arakawa  model  of  the  atmospheric  general 
circulation  is  one  of  the  more  widely  known  formulations  used  in  ex¬ 
tended  numerical  simulations  of  global  weather  and  climate.  This 
model  is  the  simplest  in  terms  of  the  physical  variables  retained, 
while  still  explicitly  treating  the  cyclone-scale  uiotlons:  it  ad¬ 
dresses  the  wind,  geopotential,  and  temperature  at  only  two  tropo¬ 
spheric  levels,  and  the  moisture  at  only  the  lower  level.  In  addi¬ 
tion,  only  the  more  important  elements  of  the  heat  and  hydrological 
balance  are  simulated  at  the  earth's  surface.  As  a  result  of  this 
simplicity,  the  two-level  Mintz-Arakawa  model  is  an  economical  for¬ 
mulation  for  producing  a  climatic  simulation,  requiring  only  26  mi¬ 
nutes  on  a  360/91  computer  to  simulate  one  day.  This  simplicity 
and  economy  would  be  irrelevant,  of  course,  were  it  not  also  a  fact 
that  this  model  has  proven  capable  of  reproducing  the  major  features 
of  the  mean  atmospheric  circulation.  For  these  reasons,  the  two- 
level  Mintz-Arakawa  model  has  been  adopted  as  the  major  simulation 
tool  of  Rand's  research  program  in  climate  dynamics,  and  a  compre¬ 
hensive  documentation  of  the  model  has  been  prepared  (Gates,  Batten, 
Kahle,  and  Nelson,  1971). 

The  atmosphere  is  driven  primarily  by  the  forcing  functions  rep¬ 
resenting  the  net  heating  and  the  net  kinetic-energy  dissipation.  To¬ 
gether  with  the  distribution  of  the  excess  evaporation  over  precipi¬ 
tation,  and  the  boundary  conditions  at  the  surface,  these  source  terms 
arc  ultimately  responsible  for  the  character  of  the  mean  circulation 
and  climate.  Since  friction,  heating,  precipitation  and  evaporation 
are  not  specified  ahead  of  time  in  this  model,  but  are  determined  by 
the  simulated  circulation  itself,  the  structure  of  the  derived  average 
forcing  fields  provides  an  important  calibration  of  the  modeling  of 
the  various  physical  processes  that  control  the  addition  and  removal 
of  energy  from  the  atmosphere.  It  is  the  purpose  of  this  Report  to 
present  the  average  source  terms,  as  simulated  by  the  two-level  model 
for  a  single  January,  and  to  compare  these  with  the  available  observa¬ 
tions  . 
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The  mean  source  terms  considered  here  are  the  averages  during  a 
30-day  integration,  corresponding  to  the  month  of  January,  during  which 
the  only  change  from  the  documented  or  reference  integration  (Gates  et 
al.,  1971)  was  a  reduction  of  10  percent  in  the  incoming  solar  radia¬ 
tion  subject  to  scattering.  This  integration  was  used  for  the  present 
analysis  for  reasons  of  convenience,  and  no  discussion  of  its  possible 
physical  significance  is  given  here.  The  average  climatic  variables 
such  as  pressure,  temperature,  cloudiness,  and  wind  simulated  during 
the  reference  or  control  experiment  itself  are  presented  elsewhere 
(Gates,  1972). 


2.  THE  FORCING  FIELDS  IN  THE  MODEL 


The  dynamics  of  the  two-level  Mintz-Arakawa  model  are  formulated 
in  terms  of  a  vertical  independent  variable  o,  which  is  a  scaled 
pressure  coordinate  given  by 

0  "  (P  -  PT>/w  (1) 

where  p  is  pressure,  pT  an  assumed  constant  pressure  at  the  tropopause 
(B  200  mb,  the  top  of  the  model  atmosphere)  and  tt  ■  p  -  p  ,  the  dif¬ 
ference  between  the  (variable)  surface  pressure  p^  and  p^,.  In  tsrms 
of  this  coordinate,  the  basic  dynamical  equations  may  be  written 

~  (irV)  +  (V  •  ff$)V  +  ^  OrdV)  +  f£  x  ttV 


+  irViJi  +  orraVtr  ■  wF  (2) 

I?  (TrcpT)  +  V  •  (*cpT$)  +  |-  (TfCpdT) 

-  ira(a  +  aV  •  Vtt  +  to)  -  ttH  (3) 

ft  +  7  *  (*<$)  +  (irqd)  -  nQ  (4) 

|I  +  7  .  (,iJ>  +  |_  (,4)  .  0  (5) 


Here  Eq.  (2)  is  the  vector  equation  of  horizontal  motion,  Eq.  (3)  is 
the  thermodynamic  energy  equation,  and  Eqs.  (4)  and  (5)  are  the  con¬ 
tinuity  equations  for  moisture  and  mass,  respectively.  Together  with 
the  equations  of  state  (o  ”  RT/p)  and  hydrostatic  balance  (3$/3o  +  tt<x 
■  0) ,  these  are  the  so-called  "primitive"  equations  of  the  model,  ex¬ 
pressed  in  'terms  of  the  horizontal  vector  velocity  $,  the  temperature  T, 
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the  specific  volume  a,  the  mixing  ratio  q,  and  the  geopotential  $, 
with  the  remaining  symbols  having  their  usual  meanings.  Of  principal 
concern  here  are  the  forcing  terms  on  the  right-hand  sides  of  Eqs. 

(2),  (3),  and  (4),  expressed  in  terms  of  £  (the  horizontal  frictional 
force  per  unit  mass),  H  (the  diabatic  heating  rate  per  unit  mass), 
and  Q  (the  rate  of  moisture  addition  per  unit  mass). 

The  boundary  conditions  accompanying  this  system  are:  a  (-  da/dt)  - 
0  at  both  0-0  and  o  ■  1,  ij>  -  gz  at  the  earth's  surface,  where  z  is 

S  g 

the  ground  elevation  (■  0  over  ocean),  and  the  assumptions  of  zero  net 
heat  flux  between  the  atmosphere  and  an  underlying  land  surface  and  of 
a  prescribed  temperature  for  an  underlying  ocean  surface.  Also  pre¬ 
scribed  in  the  model  are  the  locations  of  sea  ice  and  of  ice-covered 
land,  where  the  surface  temperature  is  not  allowed  to  rise  above  0  deg  C. 
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— — THE  global  distribution  of  the  simulated  forcing  fi ftps 

A.  THE  NET  FRICTIONAL  FORCE,  t 

In  the  application  of  Eq.  (2)  at  the  discrete  model  levels  a  -  1/4 

(level  1)  and  a  -  3/4  (level  3),  the  frictions*  force  per  unit  mass  is 
given  by 


J3  '  -fl  -  V*Vl5sl  +  «>  I8  (7) 

Here  fj  Is  an  Internal  frictional  force  at  level  1  proportional  to  the 
vertical  stress,  and  hence  proportional  to  the  vertical  shear  (3$/3z) 
at  level  2  (a  -  1/2),  with  y  an  empirical  coefficient  (-  0.44  mb  sech 
z  the  vertical  geometric  height,  and  the  factor  2g/ir  the  reciprocal  of 
the  mass  per  unit  horizontal  area  in  the  upper  and  lower  model  layers. 
The  frictional  force  ?3  at  the  lower  level  3  consists  of  a  component 
equal  and  opposite  to  F1  and  a  component  representing  the  frictional 
force  at  the  surface.  Hence  the  net  frictional  force  t  on  an  atmo¬ 
spheric  column  in  the  model  is  given  by 

'  -  ?1  +  f3  -  -VA<l5sl  +  G)  (8) 

where  CD  is  the  surface  drag  coefficient,  p 4  the  surface  air  density, 

Vs  a  measure  of  the  surface  wind  [-  0.7  V^,  where  \  (m  j  V  -  -  V  )  is 

the  extrapolated  surface  wind],  and  G  an  empirical  correction  for  lusti¬ 
ness  (■  2.0  m  sec"1). 

The  force  t  given  by  Eq.  (8)  thus  represents  the  surface  skin  fric¬ 
tion  or  alternatively,  the  force  -t  represents  the  drag  of  the  earth 
on  the  atmosphere.^  The  zonal  and  meridional  components  of  the  surface 
wind  stress  tw  h  -t  are  given,  respectively,  by 
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v-  +G)  ?  w 

tw-cdV.(^.I+g)^  <10) 

where  ufl  ■  a  cos  <p  dA/dt  is  the  zonal  surface  wind  component  and 
vfl  »  a  dcp/dt  is  the  meridional  surface  wind  component,  with  9  the 
latitude  and  A  the  longitude  on  a  spherical  earth  of  radius  a. 

In  the  model,  the  drag  coefficient  is  given  by  the  empirical 
expressions 

over  ocean: 

Cp  -  0.0025  or  (0.0010  +  0.00007|vJ  m_1  sec) 

whichever  is  smaller  (11) 

over  land  or  ice: 

C_  -  0.0020  +  0.0000012  z  m_1  (12) 

D  s 

Thus  the  surface  drag  is  velocity  dependent  over  the  oceans,  and 
varies  between  0.001  and  0.0025.  Over  land  (or  ice)  the  drag  depends 
only  on  the  surface  elevation,  and  varies  between  0.002  over  lowlands 
to  about  0.007  over  the  highest  mountains.  For  further  details,  see 
Gates  et  al.  (1971) . 

The  global  distribution  of  the  average  zonal  surface  stress 
obtained  from  the  January  simulation  is  given  in  Fig.  1,  which  clearly 
shows  the  presence  of  the  mean  surface  midlatitude  westerlies  and  the 
mean  easterlies  of  lower  latitudes  in  both  the  winter  and  summer  hemi¬ 
spheres.  These  stresses  may  be  compared  with  those  for  December- 
January-February  given  by  Hellerman  (1967) ;  the  most  systematic  dif¬ 
ferences  are  over  the  North  Pacific  ocean,  where  the  maximum  westerly 
and  easterly  stresses  are  simulated  about  15*  latitude  north  of 
their  observed  positions  and  are  approximately  50  percent  too  strong. 
Over  the  remaining  oceans,  both  the  stress  pattern  and  magnitude  are 
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in  reasonably  good  agreement  with  Hellerman' s  data.  Since  Hellerman 
used  a  quadratic  stress  formulation  similar  to  Eq.  (9)  to  determine 
the  stress  from  observed  wind-rose  data,  this  agreement  reflects  the 
general  fidelity  of  at  least  the  pattern  of  the  simulated  surface  wind 
field.  We  m»y  also  note  that  the  oceanic  drag  coefficient,  Eq.  (11), 
assumes  approximately  the  same  values  as  the  coefficient  used  by 
Hellerman. 

Over  the  continents  we  find  the  same  general  latitudinal  dependence 
of  the  simulated  zonal  stress  (Fig,  1)  as  over  the  oceans,  but  here 
the  height  dependence  of  the  drag  coefficient,  Eq.  (12),  has  produced 
local  maxima  of  the  stress  over  the  Rocky,  Andes,  and  Himalay  moun¬ 
tains,  as  well  as  systematically  increasing  the  stress  over  GreeiJand 
and  Antarctica.  The  fidelity  of  these  continental  stress  distribu¬ 
tions  over  land  is  unknown  in  the  absence  of  suitable  observations, 
although  the  height  dependence  of  in  Eq.  (12)  was  taken  from  the 
drag  coefficients  suggested  by  Cressman  (1960). 

The  corresponding  distribution  of  the  average  meridional  surface 
stress  tw  obtained  from  the  January  simulation  is  shown  in  Fig.  2. 

Here  the  agreement  over  the  oceans  with  the  data  of  Hellerman  (1967) 
is  only  fair;  the  northward  stress  over  the  northeastern  North  Pacific 
and  Atlantic  oceans  is  approximately  twice  the  strength  of  that  ob¬ 
served,  while  the  observed  broad  band  of  southward  stress  in  the  north¬ 
ern  hemisphere  subtropical  oceans  is  not  well  simulated  in  the  model 
Integration.  In  the  southern  hemisphere  the  maxima  of  northward  stress 
just  west  of  South  America,  Africa,  and  Australia  at  about  20°S  to  30°S 
correspond  well  with  observation.  Over  the  continents,  the  influence 
of  topography  is  again  noticeable  in  the  presence  of  meridional  stress 
maxima  over  the  Rockies,  Andes,  and  Himalayas,  as  well  as  over  portions 
of  Greenland  and  Antarctica.  These  average  patterns  of  meridional 
stress  reflect  the  standing  components  of  the  mean  meridional  circula¬ 
tion,  which  in  general  consists  in  the  northern  hemisphere  of  a  south¬ 
ward  flow  over  the  relatively  cold  continents  with  compensating  north¬ 
ward  flow  over  the  intermediate  warmer  oceans.  In  the  southern  hemi¬ 
sphere  the  summer  (January)  meridional  circulation  is  dominated  by  the 
subtropical  oceanic  anticyclones. 
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Although  it  is  the  stress  t  ■  F  which  serves  as  a  momentum  source 
(or  sink)  in  the  model's  equation  of  motion  (2),  the  frictional  dissi¬ 
pation  of  kinetic  energy  is  given  by  V  •  F.  This  includes  dissipation 
at  levels  1  and  3,  with  the  latter  having  a  component  dependent  upon 
the  surface  wind  V  through  Eq.  (7).  Although  it  is  thus  not  relevant  to 
the  model's  energy  budget  (to  be  considered  later),  the  surface  dissi- 

HN  “2 

pation  V  •  t  alone,  however,  has  a  global  mean  of  1.4  watts  m  in 

8  S  _2 

the  present  simulation,  which  compares  well  with  the  value  1.5  watts  m 

found  for  the  northern  hemisphere  spring  by  Ellsaesser  (1969)  using  data 

of  Kung. 

B.  THE  NET  DIABATIC  HEATING  RATE,  H 

When  the  thermodynamic  energy  equation  (3)  is  applied  at  level  1 
(o  -  1/4)  and  level  3  (o  ■  3/4),  the  diabatic  heating  rates  and 
are  given  in  terms  of  temperature  change  by 


«1-(A1  +  R2-  V^-+4I1,CM  + 

P 


AT 


1,CP 


(13) 


h3  -  (A3  +  R4  -  Kj  +  r)  I*-  +  4T  +  4T3iCP  +  AT  (14) 

p 

Here  A^  and  A^  are  the  net  short-wave  radiation  absorbed  at  levels  1 
and  3,  and  Rq,  R?,  and  R^  are  the  net  upward  long-wave  radiation  emitted 
from  the  levels  0  (o  ■  0) ,  2  (o  ■  1/2),  and  4  (o  ■  1) ,  respectively. 
Hence  A^  -f  R^  -  Rq  represents  the  net  radiation  absorbed  In  the  upper 
model  layer  (identified  with  level  1),  and  Aj  +  -  R^  is  the  net  ra¬ 

diation  absorbed  in  the  lower  layer  (identified  with  level  3).  The 
factor  2g/irc^  converts  the  heating  rates  to  the  equivalent  rate  of 
change  of  temperature.  The  term  T  represents  the  upward  flux  of  sen¬ 
sible  heat  from  the  earth's  surface,  and  is  given  by 


r 


+  G)(Tg  -  v 


(15) 
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Where  Tg  Is  the  ground  (or  ocean  surface)  temperature,  is  the  sur¬ 
face  air  temperature,  and  the  other  terms  are  as  in  Eqs.  (8)  to  (12). 

The  remaining  terms  in  Eqs.  (13)  and  (14)  are  the  rates  of  temperature 
change  due  to  the  so-called  midlevel  moist  convection  (CM),  penetrating 
moist  convection  (CP),  and  large-scale  condensation  (LS).  The  two  forms 
of  moist  convection  are  parameterized  at  both  levels  oi  the  model  in 
terms  of  the  vertical  gradients  of  the  total  heat  or  static  energy 
(cpT  +  gz  +  Lq) ,  whereas  the  large-scale  condensation  depends  upon  the 
moisture  condensed  as  a  result  of  supersaturation  at  level  3  only.  De¬ 
tails  of  these  parameterizations  are  given  elsewhere  (Gates  et  al. , 
1971). 

The  net  heating  rate  of  an  atmospheric  column  in  the  model  may 
thus  be  written 


V  -  (Hl  +  H3)  -  ai  +  a3  +  r4  -  Rq  +  r  +  ci  +  c3  +  (PREC)  Li  (16) 

where  C^  -  (AT^^  +  AT^^p)irCp/2g  is  the  upper-level  convective  heat- 
n8.  C3  -  (ATj^  +  AT3^Cp)ircp/2g  i8  the  lower-level  convective  heating, 
PREC  is  the  large-scale  condensation  rate,  and  L  is  the  latent  heat  of 
condensation. 

The  global  distribution  of  the  average  net  heating  V  obtained 
from  the  30-day  January  simulation  is  shown  in  Fig.  3.  We  note  that 
over  the  continents  of  the  northern  hemisphere  there  is  virtually  no 
average  net  heating  of  the  atmosphere.  This  may  be  influenced  by  the 
assumed  lower  boundary  condition  of  zero  net  heat  flux  at  the  earth's 
bare-land  or  ice-covered  surfaces,  i.e., 

R4  +  r  +  HE  -  Sg  -  I  .  0  (17) 

where  Hp  is  the  flux  of  latent  heat  due  to  surface  evaporation,  S  is 
the  short-wave  radiation  absorbed  at  the  land  or  ice  surface,  and8I 
is  the  upward  heat  flux  by  conduction  through  sea  ice.  This  equilib¬ 
rium  is  in  fact  used  to  determine  the  appropriate  ground  temperature, 
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T  ,  upon  which  the  terms  R, ,  r,  and  I  depend.  Thus,  over  land,  the 

O  ™ 

sensible  and  latent  heat  fluxes  are  used  in  effect  to  balance  the 
surface  heat  budget,  whereas  over  the  oceans,  where  the  surface  tem¬ 
perature  is  held  constant,  the  condition  (17)  does  not  apply.  This 
results  in  an  average  net  atmospheric  heating  where  the  prescribed 
sea-surface  temperature  is  higher  than  the  simulated  air  temperature, 
a  condition  that  occurs  over  much  of  the  world's  oceans,  as  may  be 
seen  in  Fig.  3.  This  is  particularly  marked  in  the  western  portions 
of  the  North  Atlantic  and  North  Pacific  oceans,  where  the  net  heating 
exceeds  1000  ly  day  Over  most  of  the  southern-hemisphere  oceans  the 
average  net  heating  rate  is  negative  and  of  the  order  of  -100  ly  day 
with  relatively  weak  gradients. 

The  dominant  Influence  on  the  pattern  of  net  heating  in  Fig.  3 
is  the  distribution  of  the  simulated  convective  condensation  heating, 
and  Cy  Because  these  terms  depend,  in  part,  upon  the  air-sea 
temperature  difference  and  on  the  available  moisture  in  the  air,  it 
is  not  surprising  that  the  net  heating  bears  a  resemblance  to  the 
distributions  of  both  the  sensible  heat  flux  and  the  surface  evapora¬ 
tion  (as  discussed  below) .  It  is  clear  from  these  data  that  the  oceans 
exert  an  Important  controlling  influence  on  the  model  atmosphere's  net 
heat  balance,  and  hence  on  the  simulated  mean  circulation. 

C.  THE  NET  RATE  OF  MOISTURE  ADDITION,  E  -  P 

When  the  moisture  continuity  equation  (4)  is  applied  at  level  3 
in  the  model,  the  moisture-source  term  Q  is  given  by  the  difference 
between  surface  evaporation  and  total  precipitation,  since  it  is  as¬ 
sumed  that  all  of  the  atmospheric  moisture  is  carried  at  this  single 
level .  Thus 


Q  -  Q3  -  E  -  P 

where  the  surface  evaporation  E  is  given  by 


(18) 


+  G>%  -  V 


E 


(19) 
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where  p  Is  the  density  of  water,  q.  Is  the  mixing  ratio  of  the  air 
w  ** 

juat  above  the  ground,  and  q^  is  a  mixing  ratio  for  the  ground  surface 
Itself  (which  Includes  a  measure  of  ground  wetness;  see  Gates  et  al., 
1971,  for  details).  The  total  precipitation,  P,  is  assumed  equal  to 
the  total  condensation,  and  consists  of  the  moisture  released  by  both 
the  large-scale  and  convective  processes.  Thus,  we  may  write 

P  -  PLS  +  Pc  (20) 

where  the  large-scale  precipitation  (or  condensation)  rate  is  given  by 


(PREC) 


(21) 


and  the  convective  precipitation  (or  condensation)  rate  is  given  by 


P 


C 


Cl  +  C3 


(22) 


where  PREC,  and  C3  are  as  in  Eq.  (16)  and  is  the  density  of 
water. 

In  terms  of  the  moisture  released  and  the  consequent  temperature 
change,  we  may  write 


PREC  -  Aq3jLS 


(23) 


ttL 


Ci  +  C3  "  ^3,CM  +  Aq3,CP^  2g 

ire 

“  (AT1,CM  +  AT1,CP  +  AT3,CM  +  ^.CP5  2g 


(24) 


The  moisture  temporarily  stored  in  the  air  as  cloud  is  not  considered 
in  the  present  model,  although  measures  of  ground  wetness  and  runoff 
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are  Included.  We  may  also  note  that  all  precipitation  is  assumed  to 
fall  as  rain;  snow  cover  Is  prescribed  (as  a  function  of  latitude)  only 
for  the  purpose  of  determining  the  surface  albedo. 

The  global  distribution  of  the  average  net  rate  of  moisture  addl- 
tion  Q  ■  E  -  P  obtained  from  the  30-day  simulation  is  shown  in  Fig.  4. 
This  field  reflects  the  presence  of  the  lower-latitude  precipitation 
maxima  over  the  oceans  (principally  convective  precipitation) ,  and  to 
this  extent  bears  an  inverse  resemblance  to  the  distribution  of  the 
average  net  heating  rate  in  Fig.  3.  In  the  subtropics,  the  relatively 
high  evaporation  over  the  oceans  is  the  dominant  moisture  source,  while 
in  the  higher  latitudes  the  large-scale  precipitation  serves  as  the 
principal  moisture  sink.  Over  the  North  Pacific  and  North  Atlantic 
oceans,  this  distribution  is  in  broad  agreement  with  that  derived  from 
observations  for  the  northern  winter  by  Jacobs  (1951) ,  although  his 
maxima  are  only  about  half  those  simulated  here.  Over  the  continents 
in  the  winter  hemisphere  there  is  apparently  little  net  contribution 
to  the  moisture  balance.  In  the  simulation,  as  in  nature,  the  oceans 
provide  most  of  the  moisture,  and  receive  the  bulk  of  the  precipita¬ 
tion  in  return. 


CUMULATIVE  MOISTURE  ADDITION  (MM)  ISOLINES  AT  -24.000 


The  average  net  rate  of  moisture  addition  (E  -  P)  in  a  January  simulation.  The  isolines 
are  at  2  mm  day-1  intervals,  with  the  0  isoline  dashed  (see  also  Fig.  1). 
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atnospherlc  model  should  be  corrected  before  they  are  passed  on  to  a 
coupled  ocean  circulation. 

Although  less  familiar  to  oceanographers,  the  zonal  variation 
of  the  mean  meridional  surface  win,!  stress  of  Fig.  2  could  equally 
well  drive  an  ocean  circulation,  as  the  effective  dynamical  coupling 
is  through  curl  t^.  From  Fig.  2,  however,  we  note  that  only  over  the 
North  Atlantic  would  this  stress  significantly  contribute  to  the 
(cyclonic)  ocean  circulation;  over  most  of  the  world's  oceans  the 
zonal  variation  of  (Fig.  2)  is  somewhat  smaller  than  the  meridio¬ 
nal  variation  of  tX  (Fig.  1). 

w 

B.  THE  COMPONENTS  OF  THE  NET  HEATING  RATE 

The  latitudinal  distributions  of  the  zonal  average  of  the  com¬ 
ponents  of  the  net  heating,  Eq.  (16),  are  shown  in  Figs.  6  to  9.  The 
total  latent  heat  +  (PREC)Lir/2g  released  by  the  model  's  con¬ 

densation  processes  is  shown  in  Fig.  6,  together  with  that  released 
by  the  large-scale  condensation  alone.  From  the  difference  between 
these  curves,  the  convective  processes  are  seen  to  be  the  more  impor¬ 
tant  for  total  atmospheric  heating,  -md  completely  dominate  the  heat¬ 
ing  rates  in  the  lower  latitudes.  In  the  middle  latitudes  of  both 
hemispheres,  the  convective  and  large-scale  condensation  heating  make 
approximately  equal  contributions.  North  of  about  60*N,  nearly  all 
of  the  simulated  latent  heating  is  the  result  of  large-scale  (rather 
than  convective)  condensation.  Over  the  globe,  the  average  latent 
heating  rate  is  238  ly  day"1,  of  which  201  ly  day"1  are  from  convec¬ 
tive  condensation  alone. 

The  dominance  of  large-scale  condensation  noted  above  probably 
exists  in  the  higher  latitudes  of  the  southern  hemisphere  as  well. 

The  relatively  large  amount  of  convective  latent  heating  actually  sim¬ 
ulated  between  70*S  and  90*S  (Fig.  6)  is  the  result  of  an  inadvertent 
(and  only  recently  discovered)  error  in  the  model's  treatment  of  the 
albedo  of  snow  and  ice.  This  error  has  caused  the  albedo  over  Antarc¬ 
tica  in  particular  to  be  underestimated  by  about  0.4,  giving  increased 
absorption  of  solar  radiation  at  the  surface.  This  has  in  turn  resulted 


Fig.  6  —  The  latitudinal  distribution  of  the  zonal  average  of  the  simulated  January  latent 
heating  rate  due  to  large-scale  condensation,  given  by  (PREC)Lir/2g,  and  the  total 
latent  heating  rate  due  to  both  large-scale  and  convective  condensation  (C}  +  C3) 
as  in  Eq.  (16). 
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in  anomalously  high  surface  temperatures,  from  which  the  simulation  of 
convection  has  followed  by  virtue  of  the  decreased  vertical  stability. 
The  results  of  this  albedo  error  are  also  noticeable  in  the  simulated 
fields  of  radiative  heating,  sensible  heat  flux,  precipitation,  and 
evaporation  discussed  below  (see  also  Gates,  1972). 

Figure  7  shows  the  latitudinal  distribution  of  the  net  rate  of 
short-wave  (solar)  radiation  absorption  by  water  vapor  and  clouds  in 
the  atmospheric  column;  in  Eq.  (16)  this  is  given  by  +  A3>  the  sum 
of  the  net  insolation  absorption  rate  in  the  upper  and  lower  model 
layers.  The  model's  rate  of  absorption  is  close  to,  but  systematically 
less  than,  the  observed  data  given  by  London  (1957)  for  the  northern 
hemisphere  winter.  In  the  southern  hemisphere,  the  simulated  short¬ 
wave  radiative  heating  rate  is  significantly  greater  than  that  given 
for  January  by  Sasamori,  London,  and  Hoyt  (1971).  Some  of  this  dis¬ 
crepancy,  especially  at  the  higher  southern  latitudes,  is  probably  due 
to  the  model  s  error  in  the  albedo  of  ice,  which  permits  about  twice 
as  much  solar  radiation  to  be  absorbed  at  the  surface  as  is  observed. 

Also  shown  in  Fig.  7  is  the  latitudinal  distribution  of  the  net 
rate  of  long-wave  radiative  heating,  given  by  R^  -  RQ  in  Eq.  (16),  the 
difference  between  the  long-wave  emission  from  the  surface  and  that 
from  the  top  of  the  model  atmosphere.  Since  R^  -  RQ  <  0,  the  long¬ 
wave  radiation  produces  a  net  cooling  of  the  atmospheric  column.  For 
the  January  conditions  simulated,  there  is  no  solar  radiation  absorbed 
north  of  about  70°N,  whereas  over  the  remainder  of  the  globe,  both  the 
short-  and  long-wave  radiative  heating  are  affected  by  the  simulated 
distributions  of  water  vapor  and  cloudiness .  The  agreement  with  the 
data  of  London  (1957)  and  of  Sasamori,  London,  and  Hoyt  (1971)  is  good 
north  of  30  S;  south  of  this  latitude  there  is  a  systematic  underes¬ 
timation  of  the  long-wave  cooling  rate  by  about  30  percent. 

The  maximum  cooling  rate  given  by  the  total  radiative  heating 
(Ai  +  A3  +  R^  -  Rq)  occurs  in  the  northern  subtropical  latitudes,  and 
agrees  reasonably  well  with  observation.  The  minimum  heating  simulated 
in  the  higher  latitudes  of  the  southern  hemisphere,  however,  is  an  ex¬ 
ception.  Here  the  model's  ice— albedo  error  noted  earlier  may  be  the 
cause  of  the  discrepancy  south  of  65°S.  Over  the  southern  oceans  (40*S 
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Fig.  8  —  The  latitudinal  distribution  of  the  zonal  average  of  the  simulated  January  heat 
flux  due  to  the  vertical  (turbulent)  transfer  of  sensible  heat,  as  given  by  T 
[Eq.  (15)].  The  observations  (open  circles)  are  from  Schutz  and  Gates  (1971), 
based  on  data  of  Budyko  (1963) . 


-26- 


as  given  by  V  in  Eq.  (16);  this  Is  the  sum  of  the  distributions  shown 
in  Figs.  6,  7,  and  8  (or  the  zonal  average  of  the  data  of  Fig.  3).  The 
principal  simulated  January  heating  of  the  atmosphere  is  seen  to  occur 
between  about  25®N  and  25® S,  and  is  approximately  symmetric  about  the 
equator.  Comparison  with  Fig.  6  shows  that  the  distribution  is  strongly 
influenced  by  the  latent  heat  released  through  convective  precipitation. 
In  comparison  with  the  observed  net  heating  of  the  atmosphere  given  by 
Newell  et  al.  (1970),  the  simulation  in  the  northern  hemisphere  may  be 
considered  satisfactory,  except  for  the  model's  northward  displacement, 
by  some  15°,  of  the  subtropical  zone  of  minimum  heating.  In  the  south¬ 
ern  hemisphere  the  agreement  is  much  less  satisfactory,  although  the 
simulated  and  observed  heating  curves  have  generally  similar  features. 

The  positive  total  heating  simulated  around  70° S  is  due  to  both  the  ex¬ 
cessive  short-wave  absorption  and  the  excessive  convective  latent  heat¬ 
ing  noted  earlier.  In  the  global  average,  the  computed  net  heating 
rate  is  -16  ly  day-1;  this  corresponds  to  a  cooling  of  the  entire  at¬ 
mosphere  over  the  30-day  simulation  of  about  1  deg  C,  and  is  not  in¬ 
consistent  with  the  conditions  of  this  particular  experiment  (See  Sec¬ 
tion  1) . 

C.  THE  COMPONENTS  OF  THE  MOISTURE  BALANCE 

The  latitudinal  distribution  of  the  zonal  average  of  the  simulated 
January  precipitation  is  shown  in  Fig.  10.  Here  we  see  that  the  total 
precipitation  (equal  to  the  total  condensation  in  the  model)  is  a  rea¬ 
sonably  accurate  approximation  to  the  observed  December-January-February 
precipitation  as  given  by  both  Lvovitch  and  Ovtchinnikov  (1964)  and  by 
MSller  (1951)  in  the  middle  and  high  latitudes  of  both  hemispheres.  An 
exception  occurs  over  Antarctica,  where  the  model's  simulation  is  doubt¬ 
ful  due  to  the  albedo  error  not<id  earlier).  Between  about  30°N  and  30°S, 
however,  the  observed  precipitation  is  overestimated  in  the  simulation 
by  approximately  50  percent,  with  the  tropical  rainfall  maxima  displaced 
10°  to  15°  toward  the  poles.  This  error  is  primarily  responsible  for 
the  simulated  global  average  January  precipitation  rate  of  4.11  mm  day 
being  so  far  above  the  observed  rate  of  2.94  mm  dayr^  [from  the  data  of 
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Lvovitch  and  Ovtchinnikov  (1964)]  or  2.20  mm  day  ^  [from  the  data  of 
MBller  (1951)].  Between  70#S  and  90#S,  the  anomalously  large  simulated 
precipitation  rate  is  due  to  the  convection  resulting  from  high  surface 
temperatures,  which  in  turn  are  caused  by  the  model's  inadvertent  use 
of  an  albedo  for  ice  which  is  only  about  half  the  correct  one. 

The  simulated  precipitation  in  Fig.  10  is  proportional  to  the 
total  latent  heating  in  Fig.  6,  with  a  heating  rate  of  100  ly  day"1 
corresponding  to  a  precipitation  (condensation)  rate  of  approximately 
1.72  mm  day  .  The  excessive  precipitation  modeled  in  the  tropics  is 
thus  seen  to  be  due  almost  entirely  to  convective  processes,  just  as 
was  the  total  heating  of  Fig.  9.  An  improvement  in  the  parameteriza¬ 
tion  of  convection  could  therefore  be  expected  to  aid  both  the  pre¬ 
cipitation  and  heating-rate  simulations  considerably. 

The  latitudinal  distribution  of  the  zonal  average  January  evapora¬ 
tion,  both  as  simulated  and  as  given  by  Schutz  and  Gates  (1971)  based 
on  the  data  of  Budyko  (1963),  is  given  in  Fig.  11.  As  with  the  pre¬ 
cipitation,  we  see  that  between  approximately  30°N  and  30°S  the  simu¬ 
lated  evaporation  rate  is  about  50  percent  greater  than  that  observed. 
Like  the  precipitation  error,  this  discrepancy  may  be  mainly  due  to 
convection,  with  the  surface  evaporation  attempting  to  moisten  the  at¬ 
mosphere  made  too  dry  by  excessive  rainfall.  It  is  known  that  the  sim¬ 
ulated  relative  humidity  at  level  3,  for  example,  is  only  about  half 
that  observed  at  800  mb  between  the  latitudes  30°N  and  20°S  (Gates, 
1972).  Such  a  discrepancy  would  be  sufficient  to  explain  the  evapora¬ 
tion  error,  since  the  model's  mixing  ratio  depends  primarily  upon 
the  level-3  relative  humidity  and  the  ground  wetness  [see  Gates  et  al. 
(1971)].  The  large  evaporation  rate  simulated  south  of  70°S,  like  the 
precipitation  (Fig.  10)  and  sensible  heat  flux  (Fig.  8),  is  due  to  the 
excessively  high  surface  temperatures  over  the  Antarctic  ice  (Gates, 
1972). 

The  latitudinal  distribution  of  the  zonal  average  of  the  differ¬ 
ence  E  -  P,  the  net  rate  of  moisture  addition  to  the  atmospheric  col¬ 
umn,  is  shown  in  Fig.  12  as  simulated  for  January,  together  with  two 
measures  of  the  observed  distribution.  The  simulation  may  in  general 
be  considered  a  satisfactory  approximation  to  observation.  The  sub- 
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tropics  near  30°N  are  Been  to  be  the  atmosphere's  primary  January  mois¬ 
ture  source,  with  the  region  around  30#S  a  secondary  source.  The  simu¬ 
lated  E  -  P  maximum  in  the  northern  hemisphere,  however,  is  displaced 
some  10°  north  of  the  observed,  due  to  the  excessive  precipitation 
simulated  between  10#N  and  25°N  (see  Fig.  10).  In  the  equatorial  re¬ 
gion  the  simulated  negative  E  -  P  is  only  about  half  that  observed; 
this  is  due  to  the  model's  excessive  evaporation  in  this  region,  as 
noted  earlier  (see  Fig.  11).  On  the  global  average  the  January  evap¬ 
oration  and  precipitc.'  ’on  very  nearly  balance  when  averaged  over  the 
30-day  simulation,  with  a  net  E  -  P  of  only  -0.02  mm  day”1. 
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5.  SUMMARY  AND  CONCLUSIONS 


From  the  January  simulation  summarized  here,  it  has  been  shown 
that  the  two-level  Mintz-Arakawa  model  is  capable  of  depicting  the  av¬ 
erage  forcing  or  source  fields  for  momentum,  heat,  and  moisture  with 
reasonable  accuracy.  The  simulated  surface  stress  (or  net  frictional 
force  on  an  atmospheric  column)  agrees  fairly  well  with  the  general 
pattern  of  observed  stress  derived  from  surface-wind  observations,  al¬ 
though  the  magnitude  of  the  stress  is  systematically  overestimated. 

The  simulated  net  heating  rate  of  an  atmospheric  column  shows  fair 
agreement  with  observation,  but  suffers  from  an  excessive  latent  heat¬ 
ing  rate  in  the  tropics  because  of  the  model's  apparent  overcalculation 
of  convective  condensation.  The  sensible  heat  flux  is  also  in  error  in 
the  tropics;  here  the  simulated  flux  is  in  fact  opposite  in  sign  to 
that  observed.  The  net  long-  and  short-wave  heating  rates  for  the  at¬ 
mospheric  column  agree  reasonably  well  with  observation,  but  here  too 
there  are  some  systematic  errors. 

Although  an  analysis  of  the  model's  simulation  of  the  global  Jan¬ 
uary  climate  is  given  elsewhere  (Gates,  1972),  the  present  simulation 
of  the  primary  forcing  fields  f,  H,  and  Q  of  the  system,  Eqs.  (2)  to 
(4) ,  is  sufficiently  accurate  to  lead  one  to  expect  a  reasonably  satls- 
lactory  simulation  of  the  average  climatic  elements  such  as  pressure, 
temperature,  humidity,  and  wind.  The  most  prominent  simulation  error 
which  appears  to  have  significantly  affected  the  forcing  fields  is  the 
model's  calculation  of  maximum  (convective)  precipitation  at  about  20°N 
and  20°S,  rather  than  at  the  observed  locations  near  5°N  and  5°S.  This 
error  carries  over  into  the  net  heating  rate,  and  effectively  compresses 
the  subtropical  zones  poleward  of  their  observed  positions,  especially 
in  the  northern  (winter)  hemisphere.  The  enhanced  meridional  heating 
gradient  and  the  resultant  baroclinicity  may  therefore  be  a  cause  of  the 
excessive  strength  of  the  simulated  zonal  winds,  and  of  the  associated 
excessive  strength  of  the  midlatitude  cyclones.  A  complete  analysis  of 
the  simulated  circulation's  energy  budget  and  its  partitioning  into  mean 
and  eddy  components  in  both  space  and  time  is  necessary  to  verify  this 
conjecture,  however. 
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A  second  systematic  error  noted  here  at  high  southern  latitudes 
is  the  model’s  simulation  of  too  much  absorbed  solar  radiation,  too 
much  (convective)  precipitation  and  evaporation,  and  an  upward  sen¬ 
sible  heat  flux.  These  discrepancies  are  believed  due  to  an  inadver¬ 
tent  error  in  the  albedo  of  ice,  and  are  not  considered  characteristic 
of  the  model  itself.  It  may  also  be  noted  that  the  present  results 
have  been  derived  from  a  single  30-day  January  integration,  and  not 
from  the  average  of  many  such  January  simulations.  The  year-to-year 
variations  produced  by  the  model  clearly  need  to  be  examined  before 
a  definitive  assessment  can  be  made  of  the  accuracy  of  the  simulated 
climate. 
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APPENDIX 


THE  EFFECTS  OF  SAMPLING  FREQUENCY  ON  THE  ESTIMATION  OF  THE 

MEAN  FORCING  FIELDS 


The  components  of  the  average  heating  and  moisture-addition  rates 
discussed  above  have  been  determined  from  the  30-day  net  accumulation 
of  the  various  terms,  with  the  accumulation  updated  every  time  the 
components  are  computed  in  the  program  (every  1/2  hour) .  These  fields 
have  therefore  been  determined  as  accurately  as  they  could  be  with  the 
present  model.  For  many  climatological  variables,  however,  a  less  fre¬ 
quent  sampling  is  sufficient  to  determine  a  reliable  monthly  mean.  In 
the  present  model  simulations,  a  basic  set  of  meteorological  variables 
is  normally  saved  on  a  history  tape  every  six  hours,  from  which  the 
various  components  of  the  heating,  for  example,  may  be  determined  by 
reentering  the  appropriate  portion  of  the  program  at  that  time.  As 
anticipated,  the  six-hour  sampling  frequency  was  found  to  determine 
with  sufficient  accuracy  the  simulated  January  average  pressure,  tem¬ 
perature,  and  wind  distributions  (see  Gates,  1972),  as  well  as  quanti¬ 
ties  involving  only  simple  manipulations  of  these  basic  variables.  It 
was  less  clear,  however,  that  the  net  forcing  fields  for  heat  (V)  and 
moisture  (Q)  would  be  as  accurately  portrayed  by  six-hourly  sampling. 
For  this  reason  the  averages  of  these  fields  were  determined  with  sev¬ 
eral  sampling  intervals  over  a  selected  ten-day  interval  of  the  Janu¬ 
ary  simulation,  and  the  results  are  summarized  below. 

Of  the  components  of  the  heating  rate  V,  Eq.  (16),  the  distribu¬ 
tions  of  the  net  long-wave  radiative  cooling,  R^  -  Rq,  the  sensible 
heat  flux  r,  and  the  latent  heating  due  to  large-scale  condensation 
(PREC)Lir/2g  showed  no  significant  differences  in  the  cases  of  1/2-hr 
and  6-hr  sampling  frequencies  (480  and  40  samples,  respectively,  dur¬ 
ing  the  10-day  interval) .  The  heating  rate  due  to  absorbed  solar  ra¬ 
diation,  A^  +  A^»  however,  displays  a  small  but  noticeable  change  when 
sampled  every  6  hr  rather  than  every  1/2  hr,  as  shown  in  Figs.  A1  and 
A2.  This  error  takes  the  form  of  a  spurious  enhancement  of  a  90°  long! 
tude  zonal  harmonic  tied  to  the  sun's  local  position  at  sampling  times. 


SHORT  WAVE  HEATING  ( LY )  6  -HR  SAMPLING  ISOLINES  AT  100.000  200.000 


Fig.  A1  —  The  average  rate  of  short-wave  heating  during  10  days  of  a  January  simulation 
as  given  by  sampling  every  6  hr.  The  isolines  are  at  200  ly  day”l  intervals, 
with  the  100  ly  day~^  isoline  dashed  (see  also  Fig.  1). 
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Wien  sampled  every  12  hr  (not  shown) ,  this  sampling  error  appears  as 
a  180®  longitude  harmonic  of  larger  amplitude,  and  in  the  case  of 
24-hr  sampling,  there  is  no  short-wave  absorption  at  all  over  those 
parts  of  the  earth  that  are  (always)  in  darkness  at  sampling  time. 

To  minimize  this  sampling  error,  the  product  of  the  sampling  inter¬ 
val  and  the  angular  rotation  rate  of  the  earth  (15®  hr"1)  should  not 
be  an  integral  multiple  of  the  longitudinal  grid  size  (5°).  Prefer¬ 
ably,  the  heating  rate  due  to  solar  radiation  absorption  should  be 
accumulated  each  time  it  is  computed  in  the  program.  In  Fig.  A1  the 
anomalously  large  short-wave  heating  rateB  over  Antarctica  are  due 
to  the  ice  albedo  error  noted  earlier  (see  Figs.  6  and  7).  On  the 
global  average,  the  6-hr  and  1/2-hr  sampling  give  nearly  the  same  av¬ 
erage  heating  rates  due  to  short-wave  absorption:  92.5  and  92.3 
ly  day  1,  respectively. 

Of  all  the  quantities  examined,  the  latent  heating  +  C3  in 
Eq.  (16)  due  to  the  parameterized  convective  condensation  is  the  most 
sensitive  to  sampling  frequency.  Figures  A3  and  A4  show  the  average 
convective  precipitation  rate,  Pc,  as  determined  by  both  6-hr  and  1/2-hr 
sampling  (over  the  selected  10-day  period),  to  which  the  convective 
latent  heating  rate  is  proportional  [Bee  Eq.  (22)].  In  the  case  of 
an  accurate  accumulation  of  all  convective  condensation  (1/2-hr  samp¬ 
ling)  ,  there  is  a  clear  pattern  to  the  average  convective  precipita¬ 
tion  rate  over  the  tropical  oceans;  maxima  occur  at  approximately  20®N 
and  15®S,  with  secondary  maxima  over  the  midlatitude  oceans.  In  the 
case  of  6-hr  sampling,  the  distribution  over  the  tropical  oceans  be¬ 
comes  cellular  and  confused.  This  is  due  to  the  tendency  of  the  sim¬ 
ulated  convective  precipitation  to  occur  at  isolated  points  (or  clus¬ 
ters  of  points)  for  short  periods  of  time,  usually  much  less  than  6 
hr.  Such  instantaneous  condensation  rates  should  not  be  applied  to 
periods  longer  than  those  characteristic  of  the  convection  Itself, 
which  is  here  of  the  order  of  1  hr.  Sampling  as  infrequently  as 
every  6  hr  may  also  fail  to  adequately  portray  the  diurnal  nature  of 
convection  over  the  continents  during  daytime. 

This  error  is  also  noticeable  in  terms  of  the  zonally  averaged 
distribution  shown  in  Fig.  A5.  In  the  lower  latitudes  there  is  a 


CONVECTIVE  PHEC I P 1  TAT  1  ON  6  -HH  SA^PLMG  ISOLINES  AT  0.000  10.000 


Fig.  A3  —  The  average  rate  of  convective  precipitation  during  10  days  of  a  January  simula¬ 
tion,  as  given  by  sampling  every  6  hr.  The  isolines  are  at  intervals  of  10  mm  d. 
with  the  20  mm  day“l  isoline  dashed. 
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Fig.  A5  —  The  latitudinal  distribution  of  the  zonal  average  of  the  convective  precipitation 
rate  during  10  days  of  a  January  simulation.  The  solid  curve  is  that  given  by  h - 
hr  sampling  (the  program's  frequency  of  precipitation  calculation),  and  the  dashed 
curve  that  given  by  6-hr  sampling. 
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sampling  error  of  the  order  of  25  percent  for  the  tropical  convective 
rainfall,  although  on  the  global  average,  the  6-hr  and  1/2-hr  sampling 
cases  give  nearly  the  same  average  convective  precipitation  rate  (3.52 
and  3.48  mm  day*  ,  respectively).  For  both  sampling  intervals  in 
Fig.  A2  there  is  seen  to  be  considerable  convective  precipitation  over 
Antarctica;  this  is  due  to  the  model's  inadvertent  albedo  error  for 

ice,  as  previously  discussed  in  connection  with  Figs.  6  and  10  (see 
also  Fig.  A5) . 

Although  it  shows  no  appreciable  sampling  error  as  noted  earlier, 
the  distribution  of  the  precipitation  rate  (PREC)L*/2g  due  to  large- 
scale  condensation  is  shown  in  Fig.  A6  as  determined  from  1/2-hr  ac¬ 
cumulations.  This  precipitation  is  almost  exclusively  a  middle-  and 
high-latitude  process,  and  shows  maxima  generally  in  the  western  por¬ 
tions  of  the  midlatitude  oceans.  When  added  to  the  convective  pre¬ 
cipitation  of  Fig.  A4  (1/2-hr  case),  this  distribution  gives  the  total 
precipitation  rate. 

Due  mainly  to  the  sampling  sensitivity  of  the  simulated  convec¬ 
tion  (and  the  associated  latent  heating  and  convective  precipitation), 
both  the  total  heating  rate  V  and  the  net  rate  of  moisture  addition 
E  -  P  show  a  variation  with  sampling  frequency.  If  these  rates  are 
constructed  from  data  every  6  hr  (instead  of  data  from  every  1/2  hr, 
as  in  Figs.  3  and  4) ,  a  distortion  in  low  latitudes  similar  to  that 
seen  in  Fig.  A4  occurs.  To  obtain  an  accurate  portrayal  of  the  aver¬ 
age  total-heating  and  moisture-addition  rates,  it  is  apparently  neces¬ 
sary  to  use  the  data  from  eaoh  step  of  the  calculated  convective  process 


-43- 


REFERENCES 


Budyko,  M.  I.,  Atlas  of  the  Heat  Balance  of  the  Earth,  Gidrometeorizdat, 
Moscow,  USSR,  1963,  69  pp. 

Cressman,  G.  P.,  "Improved  Terrain  Effects  in  Barotropic  Forecasts," 

Monthly  Weather  Beview,  Vol.  88,  No.  9-12,  September-December  1960, 
pp.  327-342. 

Ellsaesser,  H.  W.,  "A  Climatology  of  Epsilon  (Atmospheric  Dissipation)," 
Monthly  Weather  Beview,  Vol.  97,  No.  6,  June  1969,  pp.  415-423. 

Gates,  W.  L.,  The  January  Global  Climate  Simulated  by  the  Two-Level 
Mintz-Ardkawa  Model:  A  Comparison  with  Observation ,  R-1005-ARPA, 
The  Rand  Corporation,  Santa  Monica,  Calif.,  1972  (in  process). 

Gates,  W.  L.,  Batten,  E.  S.,  Kahle,  A.  B. ,  and  Nelson,  A.  B.,  A  Docu¬ 
mentation  of  the  Mintz-Arakawa  Two-Level  Atmospheric  General  Cir¬ 
culation  Model ,  R-877-ARPA,  The  Rand  Corporation,  Santa  Monica, 
Calif.,  1971,  408  pp. 

Hellerman,  S.,  "An  Updated  Estimate  of  the  Wind  Stress  on  the  World 
Ocean,"  Monthly  Weather  Beview,  Vol.  95,  No.  9,  September  1967, 
pp.  607-626. 

Jacobs,  W.  C.,  "The  Energy  Exchange  Between  Sea  and  Atmosphere  and  Some 
of  Its  Consequences,"  Bulletin  of  the  Saripps  Institute  of  Ocean¬ 
ography,  University  of  California,  Vol.  6,  1951,  pp.  27-122. 

London,  J.,  A  Study  of  the  Atmospheric  Heat  Balance,  Final  Report,  Con¬ 
tract  AF19(122)-165,  Department  of  Meteorology  and  Oceanography, 

New  York  University,  1957,  99  pp. 

Lvovitch,  M.  I.,  and  Ovtchinnikov,  S.  P.,  Physical-Geographical  Atlas 
of  the  World,  Academy  of  Sciences,  USSR,  and  Department  of  Geodesy 
and  Cartography,  State  Geodetic  Commission,  Moscow,  1964,  298  pp. 

MBller,  F. ,  "Vierteljahrskarten  des  Niederschlags  fllr  die  Ganze  Erde," 
Petemann's  Geographische  Mitteilungen,  Vol.  95,  No.  1,  1951, 
pp.  1-7. 

Newell,  R.  E.,  Vincent,  D.  G. ,  Dopplick,  T.  G.,  Ferruzza,  D. ,  and 

Kidson,  J.  W.,  "The  Energy  Balance  of  the  Global  Atmosphere,"  The 
Global  Circulation  of  the  Atmosphere,  Royal  Meteorological  Society, 
London,  1970,  pp.  42-90. 

Sasamorl,  T. ,  London,  J.,  and  Hoyt,  D.  V.,  Badiation  Budget  of  the 

Southern  Hemisphere,  NCAR  Manuscript  71-100,  National  Center  for 
Atmospheric  Research,  Boulder,  Colorado,  1971,  49  pp. 

Schutz,  C. ,  and  Gates,  W.  L.,  Global  Climatic  Data  for  Surface ,  800  mb, 
400  mb:  January ,  R-915-ARPA,  The  Rand  Corporation,  Santa  Monica, 
Calif.,  1971,  173  pp. 


